Confirmation of the random tiling hypothesis for a decagonal quasicrystal 
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Mechanisms that stabilize quasicrystals are much discussed but not finally resolved. We confirm 
the random tiling hypothesis and its predictions in a fully atomistic decagonal quasicrystal model by 
calculating the free energy and the phason elastic constants over a wide range of temperatures. The 
Frenkel-Ladd method is applied for the phonon part and an approach of uncorrelated phason flips 
for the configurational part. When lowering the temperature, a phase transition to an approximant 
occurs. Close to the transition temperature one of the phason elastic constants becomes soft. 

PACS numbers: 61.44.Br, 64.60.De, 63.20.Ry 



Quasicrystals are among the most complex ordered 
solid phases found in condensed matter. Their con- 
stituent atoms possess many different local environments 
and are arranged in a tiling [l[ . While perfect quasiperi- 
odic tilings like the Penrose tiling obey strict match- 
ing rules of the tiles, thermal matching rule violations 
can disorder the tiling 0, Understanding the role 
of matching rule violations for the stabilization of qua- 
sicrystals constitutes an important fundamental prob- 
lem, which so far has received attention only in idealized 
model systems 

The ensemble of matching rule violations is described 
by an order parameter, the phason strain x, which cor- 
responds to the deviation from the perfect quasipcriodic 
tiling. For small values of the phason strain, the num- 
ber of matching rule violations grows linearly, and it 
has been su gge sted that the free energy is non-analytic, 
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ground state, which is stabilized energetically, i.e. by min- 
imizing the number of energetically costly matching rule 
violations. A locked state has been confirmed with Monte 
Carlo (MC) simulations in tiling models with matching 
rule based interactions at T = in two dimensions [6[ 
and at T > in three dimensions |{|. 

In contrast, if fluctuations of phason strain are present 
in equilibrium ('unlocked'), then a quasicrystal can be 
stabilized entropically. This scenario is known as the ran- 
dom tiling hypothesis of Elser Q and Henley 0, [HJ ■ The 
free energy then follows a generalized elastic theory and 
is a quadratic form in phason strain with a minimum 
at the on average unstrained tiling, F oc \ 2 - With de- 
creasing temperature, one of the phason elastic constants 
can become soft, and a phase transition to a periodic ap- 
proximant occurs [l2| . The random tiling hypothesis has 
been investigated amply in simulations with monodis- 
perse @, @| and binary tilings [B[ and with transfer ma- 
trix methods 0, Q. Random tiling quasicrystals have 
been identified in experiment, but it is an open question 
whether they can be ground states of matter [13j. Be- 
sides quasicrystals, random tilings have been reported in 



packings of hard dimers [1J] and, recently, in molecular 
networks absorbed on graphite [15[ . 

Phonon contributions to the free energy of tilings were 
first taken into account by Koschclla et al. for the binary 
decagonal Mikulla-Roth tiling using molecular dynam- 
ics [3] and the analytic polar calculus [I?} ■ The authors 
measured the potential energy subjected to phonon strain 
and - via approximant boundary conditions - to phason 
strain at zero temperature. They showed a quadratic 
dependence of the energy on strain and determined the 
generalized elastic constants (phonon, phason, and cou- 
pling). One phason elastic constant turned out to be 
negative, rendering the quasicrystal mctastable. This 
result is in agreement with convex hull calculations of 
Lee et al. [l8| who found only crystalline ground states 
for binary tilings. The possibility of non-analytic terms 
has also been documented , but for these to exist the 
potential had to show well defined oscillations at more 
than five next neighbor distances, modeling precisely the 
matching rules. Such correlated fine tuning of a poten- 
tial, however, appears most unlikely and renders an en- 
ergetically stabilized quasicrystal improbable. 

As a rule, previous tests of the random tiling hypoth- 
esis either (i) employed discrete lattice models, (ii) dealt 
with binary tilings, which are metastable and are compli- 
cated by effects of stoichiometry and chemical potential, 
or (iii) were restricted to T = 0. Recently, an interac- 
tion pair potential with two minima, the Lennard- Jones- 
Gauss potential, was proposed [23]. At a special choice of 
positions and depths of the minima, it is possible to grow 
a stable two-dimensional monatomic decagonal random 
tiling quasicrystal from the melt in simulations. Thus a 
simple system close to reality is available, where all types 
of numerical experiments can be performed to check the 
hypotheses and paradigms of quasicrystals. Here, we 
measure the free energy in dependence of temperature 
and phason strain and derive the temperature dependent 
phason elastic constants. Our results confirm all the pre- 
dictions of the random tiling hypothesis: The free energy 
has gradient-square form, the phason elastic constants in- 
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FIG. 1. Left: A small patch of the decagonal Tubingen tri- 
angle tiling. The colors correspond to distinct vertex config- 
urations. Right: The acceptance domain in perp space is a 
regular decagon. 
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TABLE II. Phason flips change the vertex configuration of the 
Tubingen triangle tiling. The types of phason flips, difference 
between final and initial potential energy of the flipping atoms 
e, and the height of energy barrier Ea are listed. The energy 
barriers are calculated here for a single particle motion, but 
can be lowered significantly in collective motions. 
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TABLE I. The vertex configurations Vi of the Tubingen tri- 
angle tiling, numbers of neighbors in the first three shells of 
radii n = 1, r 2 = r, r 3 = 1.902, binding energies Eb and fre- 
quencies of occurrence n/N in the well-ordered state. Colors 
are identical to Fig. [1] r = (v5+ l)/2 is the golden number. 

crease monotonically with temperature, and close to the 
point where one of them vanishes, the random tiling turns 
into a rather complex but periodic approximant crystal. 

In the Lennard-Jones-Gauss system identical atoms 
interact in two dimensions by the simple double- well pair 
potential 

nr) = ^-|-eexp(-^^). (1) 

A large number of crystal structures is achieved by vary- 
ing the position r , depth e and width a of the sec- 
ond minimum. For the values r — 1.52, e — 1.8 and 
a 2 = 0.02 and at intermediate temperatures a decago- 
nal random tiling is stabilized [2(J l2l| . It is the ran- 



dom tiling variant of the Tubingen triangle tiling [22J, 
which is composed of thin and thick golden triangles and 
is related to the Penrose pattern. An alternative repre- 
sentation by polygons (regular and nonconvex decagons, 
nonagons, hexagons and pentagons) is obtained by con- 
necting nearest neighbors as shown in Fig. Q] (left). The 
Tubingen triangle tiling is relevant for modeling exper- 
imental quasicrystals including quasiperiodic colloidal 
monolayers [23j . 

By applying the cut-off radius r c = 2.0, there are nine 



tiling vertex configurations given in Tabic [I] The projec- 
tions of the vertices into the acceptance domain of per- 
pendicular (perp) space are displayed in Fig. [1] (right). 
Note that the potential energy of vertex configurations 
decreases from the center of the acceptance domain to- 
wards the outside with the particles in the center of 
decagons (dark red) being the most stable. Elementary 
excitations are local rearrangements of the tiles termed 
phason flips. The six types of phason flips are listed in 
Table |TTJ Only the vertices Vj, Vg, Vg connected to the 
boundary of the acceptance domain are able to flip. This 
is in accordance with the closure condition [24] . which 



requires that acceptance domains of random tilings are 
topologically connected when projected to perp space. 

We calculate the free energy of the decagonal tiling 
in dependence of temperature and phason strain Xij = 
djWi, where w is a two-dimensional phason displacement 
vector. Following the notation in Refs. 000 and 
in continuum harmonic approximation, the free phason 
clastic energy density per particle is 



F( X ,T)= 1 -K 2 (T) ((x$) 2 + (xg } ) 
with symmetrized phason strain modes 
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(3) 



The phason elastic constants K 2 and K\ are tempera- 
ture dependent and are weighting strain contributions of 
the Te g and T% g irreducible representations of the symme- 
try group DiQh- Both types of strain allow a symmetry 
breaking to periodic approximants of orthorhombic C 2v 
symmetry as they, when subduced, contain the trivial 
Ai-representation (r 6g |C 2 i, = A\ © A 2 = T S g\C 2v )- In 
our group theory notation we follow Ref. (25| . 
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The presence of phason strain is necessary for the 
transformation of a quasicrystal into a periodic approx- 
imant. In turn, by constructing a specific approximant 
in simulation using periodic boundary conditions, a well- 
defined average phason strain can be enforced artificially. 
We have constructed a series of orthorhombic approxi- 
mants with in the order of 20000 particles and phason 
strain |xn I , l%22 1 < 0.2, X12 = X21 = following [lf|. 
Thus, only Xeg an d Xgg were non-zero, and by varying 
them we can determine K\ and K^. 

It has been shown that typical phonon frequencies are 
at least two orders of magnitude larger than phason flip 
rates (2(| . Since the time scales of the two excitations are 
well decoupled, we can separate the free energy via F = 
-Fphon + -Fconf into a phonon part accessible to evaluation 
by molecular dynamics and a configurational part due to 
the tiling degeneracy caused by phason flips, where we 
must apply accelerated methods. 

We calculate the phonon part with a combination of 
the Frenkel-Ladd method [27[ and thermodynamic inte- 
gration as in (28|. At the low temperature To = 0.05, the 
absolute free energy is determined with Frenkel-Ladd by 
interpolating between the target structure and an Ein- 
stein crystal with known free energy. It is then extended 
to higher temperature by integrating over the internal 
energy, where we take the dependence on phason strain 
explicitly into account: 



F, 



phon 



(x,T) = E (x)+Tg(x) + I(x,T). 



(4) 



Here g(x) is an integration constant obtained by com- 
parison with Frenkel-Ladd and E{x, T) the internal en- 
ergy per particle with ground state potential energy 
-^o(x) = E(x,0)- Thermodynamic integration gives 



I(X,T) = -T AT 

J T 



, E( X ,T') - E (x) 
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For the configurational part, we apply an approxima- 
tion of uncorrelated flips. A justification for this approx- 
imation will be given below. At each value of the phason 
strain we determine the vertices with vertex configura- 
tion W, Vs, or Vg. Only these vertices are able to flip 
without introducing tiling defects and can be associated 
with a flip type Fk (see Table Hi]) . A simple Ansatz is an 
Ising-typc model with Hamiltonian 



H({s k i}) = ^2^2e k s kl 



(6) 
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using the flip energies e k (Table [TTJ) between initial state 
s k i = and final state Ski = 1, and the number of flips 
per flip type n k - The configurational free energy calcu- 
lated from the canonical partition function is 



F CO ni(x,T) = -k B T 
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FIG. 2. Contributions to the free energy for phason strain 
of symmetry Ts s - (a) Ground state potential energy Eo(x), 
(b) the Frenkel-Ladd integration constant g(x) • To, (c) the 
thermodynamic integral I(x, T) and (d) the configurational 
part -Fconf (X) T) at temperature T — 0.35. (b) and (c) are 
phason strain independent and do not contribute to the pha- 
son elastic constant, (e)-(g) The free energy at three different 
temperatures. 

Here N(x) is the number of particles in the unit cell of 
the approximant. Note that only the numbers n k depend 
on phason strain, but not the flip energies. In order to 
determine the free energy as a function of phason strain 
we measure n k for each approximant. Finally, combining 
the phonon and configuration part of the free energies 
the phason elastic constants are given by the curvature 
at zero phason strain 



K{T) 



d 2 F( X ,T) 
dx 2 



(8) 



where % = XeJ f° r A'2 and x 

Fig. [2] (a)-(d) presents the different contributions to 
the free energy at T = 0.35 for phason strain of sym- 
metry Tg g . Phason strain of symmetry Tq 9 behaves sim- 
ilarly (not shown). The part Tg(x), which arises from 
the numerically costly Frenkel-Ladd method, is large and 
shows strong fluctuations, but is practically strain inde- 
pendent. In addition, the part J(x,T), which contains 
the vibrational entropy, is also not influenced by phason 
strain. This means only the terms Eq(x) an d -Fconf (X;^) 
contribute to Eq. ([8]). They are parabolic and we have 
Ki(T) = -ftT^phon + Ki )Con f(T). The temperature depen- 
dence is contained completely in the configurational part. 

At low temperature, the configurational free energy is 
dominated by the degeneracy of the ground state (caused 
by the flips Fx, F±, F 6 , see Table HI)). The phason elas- 
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FIG. 3. Temperature dependence of the phason elastic con- 
stants K\ and K%. At low temperature, K\ is negative and 
K2 is positive. Both phason elastic constants grow linearly at 
low temperature. K\ changes sign at T\ c — 0.35 ±0.01. Error 
bars are calculated as the regression error for the parabolic 
approximation of the free energy Eq. (SJ. 

tic constants are linear in temperature, K i con f(T) oc T 
p9| . We observe that the ground state potential energy 
Eq (x) has a saddle point at \ = and the phason elastic 
constants are K\ < and K2 > (Fig. [3]), which demon- 
strates that the random tiling decagonal quasicrystal is 
not stable at T = 0. With rising temperature the con- 
figurational entropy increases both elastic constants. K\ 
changes sign at T\ c = 0.35 ± 0.01. The reversal of the 
curvature of the total free energy is shown for three tem- 
peratures in Fig. [2] (e)-(g). In agreement with this re- 
sults and close to this point, at T 2c = 0.37 ± 0.03, MC 
simulations of Ref. had already shown a phase tran- 
sition from a periodic Xi-phase to the quasicrystal. The 
Xi-phase is a (3/2, 2/1) approximant, stabilized by the 
maximization of energetically favored centered decagons, 
the lowest energy vertex (Vi) in Table fl] 

As a separate test, we calculated the free energy dif- 
ference between the Xi-phase and the unstrained decago- 
nal quasicrystal directly. The internal energy difference 
E(x,T) — Eq(x) was measured with molecular dynam- 
ics to high precision and fitted using a power series up 
to third order in T . This allowed solving the integral in 
Eq. ((5|) analytically. We find that the free energies of both 
phases cross at T-$ c = 0.39 (FigHJ), which confirms once 
more the presence of a crystal-quasicrystal transition. 

The behavior we observe resembles X-ray investiga- 
tions of icosahedral Al-Cu-Fe by Bancel et al. where 
decreasing diffraction peak intensities and a structural 
transition to a periodic phase are explained by the soften- 
ing of a phason elastic constant. In principle, a softening 
transformation can induce a single transition directly to a 
low temperature approximant or a chain of transition s by 
running through all orders of rational approximants [3l| . 
So far, the latter has not been observed in experiment. 
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FIG. 4. The free energies of the quasicrystal and the Xi- 
approximant cross at T$ c = 0.39. 



The crossing of the free energies of Xi-phase and decago- 
nal quasicrystal in FigHjpoints to a first order transition. 
A possible chain of approximants in between might have 
been suppressed because we were forced to use periodic 
boundary conditions. In the MC simulations of Ref. poj ]. 
open boundaries were used. There, indeed, the phason 
strain increases monotonically from zero to the value of 
the Xi-approximant with falling temperature, but the in- 
vestigated systems were too small to allow tracking of 
large approximants. 

It remains to justify the approximation of uncorrelated 
phason flips. We observe in the completely relaxed qua- 
sicrystal that only flips of types F\ are possible. They 
do not cost any energy, arc isolated and uncorrelated. 
Correlated flips (and eventually diffusion) arc only possi- 
ble in combination with flips of other types, which partly 
have a finite energy penalty (see Table [TTJ) and are much 
less probable (26| . Hence, the approximation of uncor- 
related flips is exact for the quasiperiodic ground state, 
and the low-temperature behavior of the phason elastic 
constants is correct to first order. In his seminal work 
on random tilings Henley commented on the existence of 
"a very close correlation between the number of flippable 
sites and the entropy per node" 



12| , a correlation that is 



analogue to our approximation of uncorrelated flips. The 
agreement of T 3c = 0.39, obtained by free energy calcu- 
lations, and the MC result T 2c = 0.37 ± 0.03 strengthens 
our confidence in the approximation of uncorrelated pha- 
son flips further. We note that the approximation does 
not affect the confirmation of the random tiling hypoth- 
esis nor any of our conclusions. 

The two-dimensional fully atomistic system studied in 
this work represents a most simple nontrivial quasicrys- 
tal. It allows to extract features which are due solely 
to quasipcriodicity. Two terms compete: (i) potential 
energy and phonon motion create a saddle point in the 
phason strain enforcing periodicity at low temperature; 
(ii) configurational entropy favors quasiperiodicity and 
dominates at high temperature. Our results agree with 
the predictions for a maximally random tiling [lij ]. al- 
though the assumption of energetic degeneracy of all 
tiling configurations is not necessary. 
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